Expression profiling of cerebrospinal fluid identifies dysregulated antiviral mechanisms in multiple sclerosis.

Despite the overwhelming evidence that multiple sclerosis is an autoimmune disease, relatively little is known about the precise nature of the immune dysregulation underlying the development of the disease. Reasoning that the CSF from patients might be enriched for cells relevant in pathogenesis, we have completed a high-resolution single-cell analysis of 96 732 CSF cells collected from 33 patients with multiple sclerosis (n = 48 675) and 48 patients with other neurological diseases (n = 48 057). Completing comprehensive cell type annotation, we identified a rare population of CD8+ T cells, characterized by the upregulation of inhibitory receptors, increased in patients with multiple sclerosis. Applying a Multi-Omics Factor Analysis to these single-cell data further revealed that activity in pathways responsible for controlling inflammatory and type 1 interferon responses are altered in multiple sclerosis in both T cells and myeloid cells. We also undertook a systematic search for expression quantitative trait loci in the CSF cells. Of particular interest were two expression quantitative trait loci in CD8+ T cells that were fine mapped to multiple sclerosis susceptibility variants in the viral control genes ZC3HAV1 (rs10271373) and IFITM2 (rs1059091). Further analysis suggests that these associations likely reflect genetic effects on RNA splicing and cell-type specific gene expression respectively. Collectively, our study suggests that alterations in viral control mechanisms might be important in the development of multiple sclerosis.

Transcript specific regulation of expression influences susceptibility to multiple sclerosis.

Genome-wide association studies (GWAS) have identified over 100 loci containing single nucleotide variants (SNVs) that influence the risk of developing multiple sclerosis (MS). Most of these loci lie in non-coding regulatory regions of the genome that are active in immune cells and are therefore thought to modify risk by altering the expression of key immune genes. To explore this hypothesis we screened genes flanking MS-associated variants for evidence of allele specific expression (ASE) by quantifying the transcription of coding variants in linkage disequilibrium with MS-associated SNVs. In total, we were able to identify and successfully analyse 200 such coding variants (from 112 genes) in both CD4+ and CD8+ T cells from 106 MS patients and 105 controls. Fifty-six of these coding variants (from 43 genes) showed statistically significant evidence of ASE in one or both cell types. In the Lck interacting transmembrane adaptor 1 gene (LIME1), for example, we were able to show that in both cell types, the MS-associated variant rs2256814 increased the expression of some transcripts while simultaneously reducing the expression of other transcripts. In CD4+ cells from an additional independent set of 96 cases and 93 controls we were able to replicate the effect of this SNV on the balance of alternate LIME1 transcripts using qPCR (p = 5 × 10-24). Our data thus indicate that some of the MS-associated SNVs identified by GWAS likely exert their effects on risk by distorting the balance of alternate transcripts rather than by changing the overall level of gene expression.

Exploring rare and low-frequency variants in the Saguenay-Lac-Saint-Jean population identified genes associated with asthma and allergy traits.

The Saguenay-Lac-Saint-Jean (SLSJ) region is located in northeastern Quebec and is known for its unique demographic history and founder effect. As founder populations are enriched with population-specific variants, we characterized the variants distribution in SLSJ and compared it with four European populations (Finnish, Sweden, United Kingdom and France), of which the Finnish population is another founder population. Targeted sequencing of the coding and non-coding immune regulatory regions of the SLSJ asthma familial cohort and the four European populations were performed. Rare and low-frequency coding and non-coding regulatory variants identified in the SLSJ population were then investigated for variant- and gene-level associations with asthma and allergy-related traits (eosinophil percentage, immunoglobulin (Ig) E levels and lung function). Our data showed that (1) rare or deleterious variants were not enriched in the two founder populations as compared with the three non-founder European populations; (2) a larger proportion of founder population-specific variants occurred with higher frequencies; and (3) low-frequency variants appeared to be more deleterious. Furthermore, a rare variant, rs1386931, located in the 3'-UTR of CXCR6 and intron of FYCO1 was found to be associated with eosinophil percentage. Gene-based analyses identified NRP2, MRPL44 and SERPINE2 to be associated with various asthma and allergy-related traits. Our study demonstrated the usefulness of using a founder population to identify new genes associated with asthma and allergy-related traits; thus better understand the genes and pathways implicated in pathophysiology.

Multiple sclerosis risk variants alter expression of co-stimulatory genes in B cells.

The increasing evidence supporting a role for B cells in the pathogenesis of multiple sclerosis prompted us to investigate the influence of known susceptibility variants on the surface expression of co-stimulatory molecules in these cells. Using flow cytometry we measured surface expression of CD40 and CD86 in B cells from 68 patients and 162 healthy controls that were genotyped for the multiple sclerosis associated single nucleotide polymorphisms (SNPs) rs4810485, which maps within the CD40 gene, and rs9282641, which maps within the CD86 gene. We found that carrying the risk allele rs4810485*T lowered the cell-surface expression of CD40 in all tested B cell subtypes (in total B cells P ≤ 5.10 × 10-5 in patients and ≤4.09 × 10-6 in controls), while carrying the risk allele rs9282641*G increased the expression of CD86, with this effect primarily seen in the naïve B cell subset (P = 0.048 in patients and 5.38 × 10-5 in controls). In concordance with these results, analysis of RNA expression demonstrated that the risk allele rs4810485*T resulted in lower total CD40 expression (P = 0.057) but with an increased proportion of alternative splice-forms leading to decoy receptors (P = 4.00 × 10-7). Finally, we also observed that the risk allele rs4810485*T was associated with decreased levels of interleukin-10 (P = 0.020), which is considered to have an immunoregulatory function downstream of CD40. Given the importance of these co-stimulatory molecules in determining the immune reaction that appears in response to antigen our data suggest that B cells might have an important antigen presentation and immunoregulatory role in the pathogenesis of multiple sclerosis.

Overexpression of the Cytokine BAFF and Autoimmunity Risk.

BACKGROUND: Genomewide association studies of autoimmune diseases have mapped hundreds of susceptibility regions in the genome. However, only for a few association signals has the causal gene been identified, and for even fewer have the causal variant and underlying mechanism been defined. Coincident associations of DNA variants affecting both the risk of autoimmune disease and quantitative immune variables provide an informative route to explore disease mechanisms and drug-targetable pathways. METHODS: Using case-control samples from Sardinia, Italy, we performed a genomewide association study in multiple sclerosis followed by TNFSF13B locus-specific association testing in systemic lupus erythematosus (SLE). Extensive phenotyping of quantitative immune variables, sequence-based fine mapping, cross-population and cross-phenotype analyses, and gene-expression studies were used to identify the causal variant and elucidate its mechanism of action. Signatures of positive selection were also investigated. RESULTS: A variant in TNFSF13B, encoding the cytokine and drug target B-cell activating factor (BAFF), was associated with multiple sclerosis as well as SLE. The disease-risk allele was also associated with up-regulated humoral immunity through increased levels of soluble BAFF, B lymphocytes, and immunoglobulins. The causal variant was identified: an insertion-deletion variant, GCTGT→A (in which A is the risk allele), yielded a shorter transcript that escaped microRNA inhibition and increased production of soluble BAFF, which in turn up-regulated humoral immunity. Population genetic signatures indicated that this autoimmunity variant has been evolutionarily advantageous, most likely by augmenting resistance to malaria. CONCLUSIONS: A TNFSF13B variant was associated with multiple sclerosis and SLE, and its effects were clarified at the population, cellular, and molecular levels. (Funded by the Italian Foundation for Multiple Sclerosis and others.).

Age Related Multiple Sclerosis Severity Score: Disability ranked by age.

BACKGROUND: The Multiple Sclerosis Severity Score (MSSS) is obtained by normalising the Expanded Disability Status Scale (EDSS) score for disease duration and has been a valuable tool in cross-sectional studies. OBJECTIVE: To assess whether use of age rather than the inherently ambiguous disease duration was a feasible approach. METHOD: We pooled disability data from three population-based cohorts and developed an Age Related Multiple Sclerosis Severity (ARMSS) score by ranking EDSS scores based on the patient's age at the time of assessment. We established the power to detect a difference between groups afforded by the ARMSS score and assessed its relative consistency over time. RESULTS: The study population included 26058 patients from Sweden ( n = 11846), Canada ( n = 6179) and the United Kingdom ( n = 8033). There was a moderate correlation between EDSS and disease duration ( r = 0.46, 95% confidence interval (CI): 0.45-0.47) and between EDSS and age ( r = 0.44, 95% CI: 0.43-0.45). The ARMSS scores showed comparable power to detect disability differences between groups to the updated and original MSSS. CONCLUSION: Since age is typically unbiased and readily obtained, and the ARMSS and MSSS were comparable, the ARMSS may provide a more versatile tool and could minimise study biases and loss of statistical power caused by inaccurate or missing onset dates.

Immunoseq: the identification of functionally relevant variants through targeted capture and sequencing of active regulatory regions in human immune cells.

BACKGROUND: The observation that the genetic variants identified in genome-wide association studies (GWAS) frequently lie in non-coding regions of the genome that contain cis-regulatory elements suggests that altered gene expression underlies the development of many complex traits. In order to efficiently make a comprehensive assessment of the impact of non-coding genetic variation in immune related diseases we emulated the whole-exome sequencing paradigm and developed a custom capture panel for the known DNase I hypersensitive site (DHS) in immune cells - "Immunoseq". RESULTS: We performed Immunoseq in 30 healthy individuals where we had existing transcriptome data from T cells. We identified a large number of novel non-coding variants in these samples. Relying on allele specific expression measurements, we also showed that our selected capture regions are enriched for functional variants that have an impact on differential allelic gene expression. The results from a replication set with 180 samples confirmed our observations. CONCLUSIONS: We show that Immunoseq is a powerful approach to detect novel rare variants in regulatory regions. We also demonstrate that these novel variants have a potential functional role in immune cells.

Increased THEMIS First Exon Usage in CD4+ T-Cells Is Associated with a Genotype that Is Protective against Multiple Sclerosis.

Multiple sclerosis is an autoimmune disease of the central nervous system. Genome wide association studies have identified over 100 common variants associated with multiple sclerosis, the majority of which implicate immunologically relevant genes, particularly those involved in T-cell development. SNP rs13204742 at the THEMIS/PTPRK locus is one such variant. Here, we have demonstrated mutually exclusive use of exon 1 and 2 amongst 16 novel THEMIS isoforms. We also show inverse correlation between THEMIS expression in human CD4+ T-cells and dosage of the multiple sclerosis risk allele at rs13204742, driven by reduced expression of exon 1- containing isoforms. In silico analysis suggests that this may be due to cell-specific, allele-dependent binding of the transcription factors FoxP3 and/or E47. Research exploring the functional implications of GWAS variants is important for gaining an understanding of disease pathogenesis, with the ultimate aim of identifying new therapeutic targets.

European multiple sclerosis risk variants in the south Asian population.

BACKGROUND: In less than a decade, genomewide association studies have identified over 100 single-nucleotide variants that are associated with increased risk of developing multiple sclerosis. However, since these studies have focused almost exclusively on European populations, it is unclear what role these variants might play in determining risk in other ethnic groups. OBJECTIVE: To assess the effects of European multiple sclerosis-associated risk variants in the south Asian population. METHODS: Using a combination of chip-based genotyping and next-generation sequencing, we have assessed 109 European-associated variants in a total of 270 cases and 555 controls from the south Asian population. RESULTS: We found that two-thirds of the tested variants (72/109) showed over representation of the European risk allele in south Asian cases (p < 0.0003). In the rest of the Immunochip array, the most associated variant was rs7318477 which maps close to TNFSF13B, the gene for the B-cell-related protein BAFF. CONCLUSION: Our data indicate substantial overlap in genetic risk architecture between Europeans and south Asians and suggest that the aetiology of the disease may be largely independent of ethnicity.

HLA associations in South Asian multiple sclerosis.

BACKGROUND: Previous efforts to identify Human Leukocyte Antigen (HLA) gene associations with multiple sclerosis (MS) in the South Asian population have been underpowered. AIM: To identify the primary HLA class II alleles associated with MS in Indians. METHODS: We typed HLA-DRB1, -DQA1 and -DQB1 in 419 patients and 451 unrelated controls by polymerase chain reaction using sequence specific oligonucleotide probes (PCR-SSOP). RESULTS: At the gene level DRB1 showed significant evidence of association (p=0.0000012), DQA1 showed only marginal evidence of association (p=0.04) and there was no evidence for association at DQB1 (p=0.26). At the DRB1 locus association is confirmed with the *15:01 (p=0.00002) and the *03 (p=0.00005) alleles. CONCLUSION: Our study confirms that the risk effects attributable to the HLA- DRB1*15:01and DRB1*03 alleles seen in Europeans are also seen in Indians. The absence of any evidence of association with DQB1 alleles reflects the lower linkage disequilibrium between DQB1 alleles and DRB1 risk alleles present in this population, and illustrates the potential value of fine mapping signals of association in different ethnic groups.

Class II HLA interactions modulate genetic risk for multiple sclerosis.

Association studies have greatly refined the understanding of how variation within the human leukocyte antigen (HLA) genes influences risk of multiple sclerosis. However, the extent to which major effects are modulated by interactions is poorly characterized. We analyzed high-density SNP data on 17,465 cases and 30,385 controls from 11 cohorts of European ancestry, in combination with imputation of classical HLA alleles, to build a high-resolution map of HLA genetic risk and assess the evidence for interactions involving classical HLA alleles. Among new and previously identified class II risk alleles (HLA-DRB1*15:01, HLA-DRB1*13:03, HLA-DRB1*03:01, HLA-DRB1*08:01 and HLA-DQB1*03:02) and class I protective alleles (HLA-A*02:01, HLA-B*44:02, HLA-B*38:01 and HLA-B*55:01), we find evidence for two interactions involving pairs of class II alleles: HLA-DQA1*01:01-HLA-DRB1*15:01 and HLA-DQB1*03:01-HLA-DQB1*03:02. We find no evidence for interactions between classical HLA alleles and non-HLA risk-associated variants and estimate a minimal effect of polygenic epistasis in modulating major risk alleles.

Dissection of a Complex Disease Susceptibility Region Using a Bayesian Stochastic Search Approach to Fine Mapping.

Identification of candidate causal variants in regions associated with risk of common diseases is complicated by linkage disequilibrium (LD) and multiple association signals. Nonetheless, accurate maps of these variants are needed, both to fully exploit detailed cell specific chromatin annotation data to highlight disease causal mechanisms and cells, and for design of the functional studies that will ultimately be required to confirm causal mechanisms. We adapted a Bayesian evolutionary stochastic search algorithm to the fine mapping problem, and demonstrated its improved performance over conventional stepwise and regularised regression through simulation studies. We then applied it to fine map the established multiple sclerosis (MS) and type 1 diabetes (T1D) associations in the IL-2RA (CD25) gene region. For T1D, both stepwise and stochastic search approaches identified four T1D association signals, with the major effect tagged by the single nucleotide polymorphism, rs12722496. In contrast, for MS, the stochastic search found two distinct competing models: a single candidate causal variant, tagged by rs2104286 and reported previously using stepwise analysis; and a more complex model with two association signals, one of which was tagged by the major T1D associated rs12722496 and the other by rs56382813. There is low to moderate LD between rs2104286 and both rs12722496 and rs56382813 (r2 ≃ 0:3) and our two SNP model could not be recovered through a forward stepwise search after conditioning on rs2104286. Both signals in the two variant model for MS affect CD25 expression on distinct subpopulations of CD4+ T cells, which are key cells in the autoimmune process. The results support a shared causal variant for T1D and MS. Our study illustrates the benefit of using a purposely designed model search strategy for fine mapping and the advantage of combining disease and protein expression data.

Statistical colocalization of genetic risk variants for related autoimmune diseases in the context of common controls.

Determining whether potential causal variants for related diseases are shared can identify overlapping etiologies of multifactorial disorders. Colocalization methods disentangle shared and distinct causal variants. However, existing approaches require independent data sets. Here we extend two colocalization methods to allow for the shared-control design commonly used in comparison of genome-wide association study results across diseases. Our analysis of four autoimmune diseases--type 1 diabetes (T1D), rheumatoid arthritis, celiac disease and multiple sclerosis--identified 90 regions that were associated with at least one disease, 33 (37%) of which were associated with 2 or more disorders. Nevertheless, for 14 of these 33 shared regions, there was evidence that the causal variants differed. We identified new disease associations in 11 regions previously associated with one or more of the other 3 disorders. Four of eight T1D-specific regions contained known type 2 diabetes (T2D) candidate genes (COBL, GLIS3, RNLS and BCAR1), suggesting a shared cellular etiology.

Long-range epigenetic regulation is conferred by genetic variation located at thousands of independent loci.

The interplay between genetic and epigenetic variation is only partially understood. One form of epigenetic variation is methylation at CpG sites, which can be measured as methylation quantitative trait loci (meQTL). Here we report that in a panel of lymphocytes from 1,748 individuals, methylation levels at 1,919 CpG sites are correlated with at least one distal (trans) single-nucleotide polymorphism (SNP) (P<3.2 × 10(-13); FDR<5%). These trans-meQTLs include 1,657 SNP-CpG pairs from different chromosomes and 262 pairs from the same chromosome that are >1 Mb apart. Over 90% of these pairs are replicated (FDR<5%) in at least one of two independent data sets. Genomic loci harbouring trans-meQTLs are significantly enriched (P<0.001) for long non-coding transcripts (2.2-fold), known epigenetic regulators (2.3-fold), piwi-interacting RNA clusters (3.6-fold) and curated transcription factors (4.1-fold), including zinc-finger proteins (8.75-fold). Long-range epigenetic networks uncovered by this approach may be relevant to normal and disease states.

General and specific effects of early-life psychosocial adversities on adolescent grey matter volume.

Exposure to childhood adversities (CA) is associated with subsequent alterations in regional brain grey matter volume (GMV). Prior studies have focused mainly on severe neglect and maltreatment. The aim of this study was to determine in currently healthy adolescents if exposure to more common forms of CA results in reduced GMV. Effects on brain structure were investigated using voxel-based morphometry in a cross-sectional study of youth recruited from a population-based longitudinal cohort. 58 participants (mean age = 18.4) with (n = 27) or without (n = 31) CA exposure measured retrospectively from maternal interview were included in the study. Measures of recent negative life events (RNLE) recorded at 14 and 17 years, current depressive symptoms, gender, participant/parental psychiatric history, current family functioning perception and 5-HTTLPR genotype were covariates in analyses. A multivariate analysis of adversities demonstrated a general association with a widespread distributed neural network consisting of cortical midline, lateral frontal, temporal, limbic, and cerebellar regions. Univariate analyses showed more specific associations between adversity measures and regional GMV: CA specifically demonstrated reduced vermis GMV and past psychiatric history with reduced medial temporal lobe volume. In contrast RNLE aged 14 was associated with increased lateral cerebellar and anterior cingulate GMV. We conclude that exposure to moderate levels of childhood adversities occurring during childhood and early adolescence exerts effects on the developing adolescent brain. Reducing exposure to adverse social environments during early life may optimize typical brain development and reduce subsequent mental health risks in adult life.

Multiple sclerosis genetics.

Genome-wide association studies have revolutionised the genetic analysis of multiple sclerosis. Through international collaborative efforts involving tens of thousands of cases and controls, more than 100 associated common variants have now been identified. These variants consistently implicate genes associated with immunological processes, overwhelmingly lie in regulatory rather than coding regions, and are frequently associated with other autoimmune diseases. The functional implications of these associated variants are mostly unknown; however, early work has shown that several variants have effects on splicing that result in meaningful changes in the balance between different isoforms in relevant tissues. Including the well established risk attributable to variants in genes encoding human leucocyte antigens, only about a quarter of reported heritability can now be accounted for, suggesting that a substantial potential for further discovery remains.

Analysis of immune-related loci identifies 48 new susceptibility variants for multiple sclerosis.

Using the ImmunoChip custom genotyping array, we analyzed 14,498 subjects with multiple sclerosis and 24,091 healthy controls for 161,311 autosomal variants and identified 135 potentially associated regions (P < 1.0 × 10(-4)). In a replication phase, we combined these data with previous genome-wide association study (GWAS) data from an independent 14,802 subjects with multiple sclerosis and 26,703 healthy controls. In these 80,094 individuals of European ancestry, we identified 48 new susceptibility variants (P < 5.0 × 10(-8)), 3 of which we found after conditioning on previously identified variants. Thus, there are now 110 established multiple sclerosis risk variants at 103 discrete loci outside of the major histocompatibility complex. With high-resolution Bayesian fine mapping, we identified five regions where one variant accounted for more than 50% of the posterior probability of association. This study enhances the catalog of multiple sclerosis risk variants and illustrates the value of fine mapping in the resolution of GWAS signals.

MANBA, CXCR5, SOX8, RPS6KB1 and ZBTB46 are genetic risk loci for multiple sclerosis.

A recent genome-wide association study reported five loci for which there was strong, but sub-genome-wide significant evidence for association with multiple sclerosis risk. The aim of this study was to evaluate the role of these potential risk loci in a large and independent data set of ≈ 20,000 subjects. We tested five single nucleotide polymorphisms rs228614 (MANBA), rs630923 (CXCR5), rs2744148 (SOX8), rs180515 (RPS6KB1), and rs6062314 (ZBTB46) for association with multiple sclerosis risk in a total of 8499 cases with multiple sclerosis, 8765 unrelated control subjects and 958 trios of European descent. In addition, we assessed the overall evidence for association by combining these newly generated data with the results from the original genome-wide association study by meta-analysis. All five tested single nucleotide polymorphisms showed consistent and statistically significant evidence for association with multiple sclerosis in our validation data sets (rs228614: odds ratio = 0.91, P = 2.4 × 10(-6); rs630923: odds ratio = 0.89, P = 1.2 × 10(-4); rs2744148: odds ratio = 1.14, P = 1.8 × 10(-6); rs180515: odds ratio = 1.12, P = 5.2 × 10(-7); rs6062314: odds ratio = 0.90, P = 4.3 × 10(-3)). Combining our data with results from the previous genome-wide association study by meta-analysis, the evidence for association was strengthened further, surpassing the threshold for genome-wide significance (P < 5 × 10(-8)) in each case. Our study provides compelling evidence that these five loci are genuine multiple sclerosis susceptibility loci. These results may eventually lead to a better understanding of the underlying disease pathophysiology.

Confounding underlies the apparent month of birth effect in multiple sclerosis.

OBJECTIVE: Several groups have reported apparent association between month of birth and multiple sclerosis. We sought to test the extent to which such studies might be confounded by extraneous variables such as year and place of birth. METHODS: Using national birth statistics from 2 continents, we assessed the evidence for seasonal variations in birth rate and tested the extent to which these are subject to regional and temporal variation. We then established the age and regional origin distribution for a typical multiple sclerosis case collection and determined the false-positive rate expected when comparing such a collection with birth rates estimated by averaging population-specific national statistics. RESULTS: We confirm that seasonality in birth rate is ubiquitous and subject to highly significant regional and temporal variations. In the context of this variation we show that birth rates observed in typical case collections are highly likely to deviate significantly from those obtained by the simple unweighted averaging of national statistics. The significant correlations between birth rates and both place (latitude) and time (year of birth) that characterize the general population indicate that the apparent seasonal patterns for month of birth suggested to be specific for multiple sclerosis (increased in the spring and reduced in the winter) are expected by chance alone. INTERPRETATION: In the absence of adequate control for confounding factors, such as year and place of birth, our analyses indicate that the previous claims for association of multiple sclerosis with month of birth are probably false positives.